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ABSTRACT 

We use our very deep Keck Deep Fields UnGTZI catalog of z~4, 3, and 2 UV-selected star-forming 
galaxies to study the evolution of the rest-frame 1700A luminosity density at high redshift — a study 
that is motivated by our finding of luminosity-dependent evolution of the galaxy luminosity function at 
high redshift. Ours is the most robust UV luminosity density measurement to date at these redshifts 
as it uses a well-tested object selection technique, several independent sightlines, and probes deep into 
the galaxy luminosity function. The ability to reliably constrain the contribution of faint galaxies 
is critical and our data do so as they reach deep into the galaxy population, to even at 

z^A and deeper still at lower redshifts (where Af£^Q=— 21.0 is the Schechter function's characteristic 
magnitude for Lyman break galaxies (LBGs) and i^SG corresponding luminosity). We find that 

the luminosity density at high redshift is dominated by the hitherto poorly studied galaxies fainter 
than i^BG' Slid, indeed, the the bulk of the UV light in the high-z Universe comes from galaxies in 
the rather narrow luminosity range L=0.1-1L2bg- these faint galaxies that govern the behavior 
of the total UV luminosity density. Overall, there is a gradual rise in luminosity density starting at 
2~4 or earlier (we find twice as much UV light at z^-^S as at z^A), followed by a shallow peak or 
a plateau within 2:^3-1, and then followed by the well- know plunge at lower redshifts. Within this 
total picture, luminosity density in sub-L^^^ galaxies evolves more rapidly at high redshift, z'>2, 
than that in more luminous objects. However, this is reversed at lower redshifts, z<l, a reversal that 
is reminiscent of galaxy downsizing, albeit here thus far seen only in galaxy luminosity and not yet in 
galaxy mass. Within the context of the models commonly used in the observational literature, there 
seemingly aren't enough faint or bright LBGs to maintain ionization of intergalactic gas even as late 
as 2:~4. This is particularly true at earlier epochs and even more so if the faint-end evolutionary 
trends we observe at Z'^S and 4 continue to higher redshifts. Apparently the Universe must be easier 
to reionize than some recent studies have assumed. Nevertheless, snh-L'^gQ galaxies do dominate 
the total UV luminosity density at z>2 and this dominance further highlights the need for follow-up 
studies that will teach us more about these very numerous but thus far largely unexplored systems. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
starburst 



1. INTRODUCTION 

The global, volume-averaged UV luminosity density of 
the Universe, and the global star formation rate density 
that is often derived from it, contain important informa- 
tion regarding the formation and evolution of galaxies, 
the populating of the Universe with stars, and the pro- 
duction of heavy elements in the cosmological context. 

The first measurement of the time evolution of the cos- 
mic UV luminosity density was done by Lilly et al. (1996) 
who found a remarkable drop in the luminosity density 
over the redshift interval from z~l to z~0 which, with 
the help of stellar population synthesis models, they in- 
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terpreted as an order-of-magnitude decline in cosmic star 
formation rate since 2~1. The public release of Hubble 
Deep Field data soon thereafter (HDF; Williams et al. 
1996) made it possible to extend such studies to higher 
redshifts, z~l-4, and in nearly concurrent work (the two 
papers were published only a month apart) Madau et al. 
(1996) placed lower limits^ on the star formation rate 
density at z^3 and 4, while Sawicki, Lin, & Yee (1997) 
estimated the UV luminosity density in M[/v(Ai?)> — 15 
galaxies over z~ 1.5-3. 5, and concluded that the onset of 
star formation in the Universe occurred at z>3.5. 

Much subsequent work has focused on either extending 
the analysis to higher redshifts (e.g., Lehnert & Bremer 
2003; Bunker et al. 2004; Yan et al. 2004; Bouwens et 
al. 2006) in attempts to identify the galaxy population 
responsible for cosmic reionization, or on assessing the 

The Madau et al. (1996) analysis included neither incomplete- 
ness corrections nor a uniform limiting magnitude. 
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corrections for dust either by measuring attenuation in 
optically-selected galaxies (e.g., Sawicki & Yee 1998; Pa- 
povich et al. 2001; Shapley et al. 2001; Reddy & Stei- 
del 2004; Iwata, Inoue, & Burgarella 2005a), tracking 
down the contribution of rare but vigorously star- forming 
sub-millimeter-selected galaxies (e.g., Small et al. 1997; 
Barger et al. 1998; Bales et al. 1999; Chapman et al. 
2005; Sawicki & Webb 2005), or attempting to relate 
these two star- forming populations (e.g.. Chapman et al. 
2000; Sawicki 2001; Baker et al. 2001, Webb et al. 2003). 

However, the early HDF-based results (e.g., Madau et 
al. 1996, 1998; Sawicki et al. 1997; Connolly et al. 1997) 
on which the picture at z~l-4 was based suffered from 
a serious drawback: they were based on a single, tiny 
field and thus subject to both Poisson noise and poten- 
tial systematic problems related to large-scale structure 
effects (often referred to as "cosmic variance"). Recog- 
nizing this as a serious problem, Steidel et al. (1999) 
used their large, spectroscopically-calibrated and multi- 
field Lyman break galaxy samples to constrain the UV 
luminosity function and cosmic star formation rate den- 
sity at 2;~3 and 4. Their conclusion was that there was 
no evidence for evolution in either the galaxy luminos- 
ity function or cosmic SFR density from to z~4. 
However, the Steidel et al. ground-based surveys contain 
only galaxies brighter than AI*+1 at z~3 and Af*at z~4. 
Consequently, while the bright end of the galaxy popula- 
tion in their study is reasonably well constrained, these 
authors, too, had to rely on HDF data to measure the 
LP below ^L*. Specifically, at they used the HDF 
LBG counts to constrain the LP's faint end, reporting it 
to be a = —1.6; however, lacking sufRcient statistics for 
faint z^A galaxies in the HDP they were then forced to 
assume the same a = —1.6 faint end slope for their z~4 
measurement. Thus, the Steidel et al. (1999) z~3 lumi- 
nosity density measurement is in fact highly reliant on 
the single tiny HDP, while their 2^4 luminosity density 
measurement, and the conclusion that the SPR density 
does not evolve between 2^4 and 2^3 is based on the 
assumption of a non-evolving faint end of the LP. 

Meanwhile, the contribution of faint, sub-_L*galaxies 
to the cosmic luminosity density and star formation rate 
cannot be stressed strongly enough. As Pig.nillustrates, 
for most reasonable luminosity functions it is exactly the 
faint, hitherto poorly constrained sub-L* galaxies that 
produce the bulk of the total light. The Steidel et al. 
(1999, 2003) ground-based surveys are so shallow that 
they capture only about a third of the total light at 2~3 
and less than ~20% at 2~4 and similar limitations apply 
to much of the recent work at higher redshifts. The infer- 
ence about the remaining yet dominant bulk of the UV 
luminosity density is based on poorly-constrained HDP- 
based extrapolations to fainter luminosities. 

Recognizing the limitations of the Steidel et al. (1999, 
2003, 2004) surveys, we have recently undertaken an 
imaging survey specifically targeting intrinsically faint, 
sub-L* Lyman break galaxies. Our survey, the Keck 
Deep Pields (KDP; Sawicki & Thompson 2005, 2006) se- 
lects 2;~4, 3, and 2 galaxies using the very same UnGTZI 
filter set and color-color selection criteria as those used 
by the Steidel team and, indeed, our photometry is even 
further calibrated onto that of Steidel et al. since some 
of our KDP fields are co-located with theirs. This ap- 
proach gives us a uniquely well-defined, characterized. 



and robust sample, but one that probes 1.5 magnitudes 
deeper than the Steidel et al. work and so reaches signif- 
icantly below L*: even at 2~4 we reach M*+2 and so 
can constrain the shape of the all-important faint end of 
the luminosity function very well. 

Using these very deep KDP data, we found a shal- 
lower LP faint-end slope at 2~3 than did Steidel et al. 
(1999): a = -1.43 vs. their HDP-based a = -1.6; and 
an even shallower one at z~4: a = —1.26 vs. their as- 
sumed a = —1.6 (Sawicki & Thompson 2006). Because 
of our Steidel-like sample selection combined with an un- 
precedented wide range of tests that we carried out on 
our sample, our LP measurements are uniquely robust 
against systematic effects. Details of the KDP obser- 
vations and LP measurements are given in Sawicki & 
Thompson (2005) and Sawicki & Thompson (2006), re- 
spectively, and are briefly revied in § |21 The purpose of 
the present paper is to explore how these new, more ac- 
curate, shallower LPs impact our view of the cosmic UV 
luminosity density and star formation history of the Uni- 
verse, by taking proper account of the evolving number 
density of faint, sub-L*galaxies. 

As in all the papers in the KDP series, in the present 
paper we use the AB flux normalization (Oke, 1974) and 
adopt f7A/=0.3, f^A=0.7, and i?o=70 km s-^Mpc-^ We 
also deflne the characteristic magnitude 7\/£gQ=— 21.0, 
and a corresponding characteristic luminosity L*j^bq- 
These paramaters correspond very closely to the value 
of Al* at 2~3 and 2^4 as measured by us in the KDP 
(Sawicki & Thompson 2006) and at 2'--^3 by Steidel et al. 
(1999). 

2. THE KDF GALAXY SAMPLE AND LUMINOSITY 
FUNCTIONS 

In this section we briefly summarize some of the most 
important features of the KDP data on which we base 
our analysis. A detailed description of the Keck Deep 
Pield observations, data reductions, and high-z galaxy 
selection can be found in Sawicki & Thompson (2005), 
while our z~4, 3, and 2 luminosity function measure- 
ments are presented and discussed in detail in in Sawicki 
& Thompson (2006). 

2.1. The KDF survey 

Our analysis is based on results from our very deep 
multicolor imaging survey carried out using a total of 
71 hours of integration on the Keck I telescope. These 
Keck Deep Pields use the same UnGTZI fllter set and 
color-color selection techniques as are used for brighter 
LBGs in the work of Steidel et al. (1999, 2003, 2004). In 
contrast to the Steidel et al. work, however, the KDPs 
reach TZum—'^.T; this is 1.5 magnitudes deeper than Stei- 
del et al. surveys and significantly below L* at 2=2-4: 
even at 2~4 we reach two magnitudes fainter than M* . 
Purthermore, the KDPs have an area of 169 arcmin^ di- 
vided into 3 spatially-independent patches on the sky — 
an arrangement that allows us to monitor the effects of 
cosmic variance on our results. 

An important feature of the KDP is that our use of 
the UnGTZI fllter set lets us select high-2 galaxies using 
the color-color selection criteria deflned and spectroscop- 
ically tested by Steidel et al. (1999, 2003, 2004). We can 
thus confidently select sub-L* star-forming galaxies at 
high redshift without the need for spectroscopic charac- 
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terization of the sample — characterization that would 
be extremely expensive at the faint magnitudes that in- 
terest us. Moreover, this commonality with the Steidel 
et al. surveys has allowed us to reliably combine our data 
with theirs, thereby for the first time consistently span- 
ning a large range in galaxy luminosity at high redshift. 

To 7^^im=27, the KDF contain 427, 1481, 2417, and 
2043 J7„G7?./-selected star-forming galaxies at z~4, z^3, 
z^2.2, and z^l.7, respectively, selected using the Stei- 
del et al. (1999, 2003, 2004) color-color selection criteria. 
At our completeness limit, the KDF data probe galax- 
ies with UV luminosities that correspond to star forma- 
tion rates (uncorrected for dust obscuration) of only 1-2 
Me/yr. 

2.2. The z^2, 3, and 4 luminosity functions 

Following Steidel et al. (1999), our LF calculation uses 
the effective volume technique, computing Veff via re- 
covery tests of artificial high-z galaxies implanted into 
the images. To increase the luminosity range over which 
we can robustly study the high-z LFs, at the bright end 
of the z~3 and z~4 samples we augmented our KDF data 
with the results of Steidel et al. (1999). Because of the 
identical filter set, object selection, and LF measurement 
technique, we can do this with little fear of introducing 
systematic offsets as can happen when other deep but rel- 
atively small-area surveys are augmented with shallower, 
wide-field data. 

In Fig. |2]we reproduce the z~4, 3, and 2.2 luminosity 
functions from Sawicki & Thompson (2006). In that pa- 
per we have examined a wide range of potential sources 
of systematic problems, including effects due to modeling 
of survey volume, differential sample selection, field-to- 
field variance, uncertainties in /c-corrections, and so on. 
On the basis of these tests we have a very good under- 
standing of the reliability of our LF results. In summary, 
wc found that at z~4 and z~3, our LF measurements are 
immune to a range of potential systematic effects, while 
our Z'--^2.2 and z~1.7 LFs are more uncertain. The z~1.7 
LF is likely quite strongly affected by systematics and so 
we do not include it in the analysis presented here. The 
Z'--^2.2 LF measurement is possibly also affected by sys- 
tematics but to a much smaller degree than the z~1.7 
one: we can confidently state that the number density 
of Z'^2.2 galaxies is not lower than that shown in Fig. El 
and that if it is higher, it is higher by no more than a fac- 
tor of ~2. The potential systematic error at z~2.2 is due 
to our uncertainty regarding the survey volume which it- 
self stems from our uncertainty in the model colors of 
z~2.2 galaxies. The LF shown in Fig.|5]is based on our 
best and most realistic assumptions, namely that there is 
moderately strong dust obscuration of E{B — y)=0.15. 
The potentially twice-higher normalization would ensue 
in the extreme and rather unrealistic case of no dust (see 
Sawicki & Thompson 2006 for details) — a scenario that 
appears highly unlikely in view of recent spectral energy 
distribution studies (Shapley et al. 2005; Sawicki et al. 
2006). It must be stressed again, however, that the z~3 
and z~4 LFs, are very robust against such and other 
systematics. 

The uncertainties in Fig. 13 include both \/N statisti- 
cal uncertainty in galaxy numbers and a measure of the 
effects of large scale structure ("cosmic variance") ob- 
tained from a bootstrap resampling analysis of the mul- 



tiple fields in the KDF. Both these sources of uncertainty 
are propagated into our error estimates for the luminos- 
ity density. We also note that at both z~4 and ^^^3, 
we find M*fii —21.0, in good agreement with the z~3 
Steidel et al. 1999 value converted into our cosmology. 
For convenience, we have thus defined a reference char- 
acteristic magnitude M£^q= —21.0 and the equivalent 
characteristic luminosity i^sG' these values correspond 
to a galaxy with a (dust-free) star formation rate of 15 
Mg/yr (e.g., Kennicut 1998). 

In contrast to the Steidel et al. (1999) results, how- 
ever, we find a somewhat shallower faint-end slope a. At 
z~3 we find a — —1.43, compared to their a = —1.6, 
and at z~4 we find a=— 1.26. This evolving faint-end 
slope marks one of the most striking findings of the KDF 
to date: a differential, luminosity-dependent evolution 
of the Lyman Break Galaxy luminosity function. Sim- 
ply put, while the number density of luminous galaxies 
{L>L*) remains fixed, the number density of sub-L* ob- 
jects more than doubles from z~4 to Z'--^3. Our numerous 
tests show that this result is highly robust and insensitive 
to systematic biases (see Sawicki & Thompson 2006). In 
the present paper we turn to investigate how this differ- 
entially evolving galaxy population, refiected in an evolv- 
ing luminosity function, impacts our view of the cosmic 
luminosity density and star formation rate density. 

3. THE LUMINOSITY DENSITY 

The luminosity density is the sum of all the light 
at a given wavelength, or range of wavelengths, from 
all galaxies in a unit volume. It can be computed as 
the luminosity-weighted integral of the galaxy luminos- 
ity function 4>{L): 

/•OO 

PL = / L<t>{L)dL, (1) 

which, for a luminosity function of the Schechter (1976) 
form, (l){L)dL = (j)*{L/ L*)°'e-^l^' dL/ L, can be ex- 
pressed in terms of the gamma function, r(a;) = 
Jg°° e~^t^~^dt, as 

PL = rL*Tia + 2). (2) 

Because of the difficulty of constraining the number 
density of faint galaxies at high redshift, most studies 
focus on the partial, or incomplete, luminosity density 
— i.e., the luminosity density due to objects brighter 
than some limiting luminosity Lum- We, too, will do 
so here at first by focusing on galaxies brighter than 
O.lL^sQ (§ 13. 1() before presenting the total luminosity 
density (§ 13.2(1 and then examining the contributions to 
that total that are made by galaxies with luminosities in 
various luminosity bins 13.3(1 . 

3.1. Luminosity in galaxies brighter than O.IL^^q 

Given a galaxy luminosity function ipiL), the luminos- 
ity density due to galaxies brighter than Lum is 

/>oo 

Pl{L > L™„) = L4>{L)dL. (3) 

For the Schechter LF parametrization this can be ex- 
panded as 

Pl{L > Lrmn) = ^1^*1^0^ + 2) - 7(0 + 2,L„„„/L*)], 

(4) 
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where j{x,l) = Jl^e~*t^~^dt is the lower incomplete 
gamma function (e.g., Arfken 1985, Press et al. 1986). 

We adopt Lum—O.lLlgQ, which corresponds to 
-^i700=~18.5 or, in the absence of dust, a galaxy with a 
star formation rate of ~1.5 Mq/jt (Kennicutt 1998). As 
Fig- n shows, this choice of Lht™ captures approximately 
75% of the UV light at z^A and 3 for the luminosity 
functions we measure in the KDF. 

We show the pL{L>0.lLlg(^) KDF values with filled 
circles in Fig. O together with GALEX values at lower 
redshift (filled squares) and Iwata et al. (2006a, 2006b) 
values at z^5 (filled triangle). Tables ^ and |2] list the 
values of luminosity density for the KDF and for other 
surveys. The KDF values were computed by applying 
Eq. O to the Schechter function parameters of Sawicki 
& Thompson (2006). The uncertainties associated with 
our KDF points in Fig. O (and in all subsequent fig- 
ures) include the contributions of both Poisson statistics 
and cosmic variance, the latter estimated through boot- 
strap resampling of the five KDF fields. Also shown in 
Fig.Oare the z<l GALEX UV luminosity density values 
of Schiminovich et al. (2005) which we have translated 
from their 1500A to 1700A and adjusted to represent 
only galaxies brighter than O.li^^g by scaling their to- 
tal luminosity densities using our Equations |31 and ^ to- 
gether with the GALEX LF Schechter parameters given 
by Arnouts et al. (2005) and Wyder et al. (2005). The 
z~5 point is based on new work by Iwata et al. (2006a, 
2006b), who select z^b LBGs from Subaru Suprimecam 
data (the median redshift of their sample is z=4.8, but 
for simplicity we refer to it here as z~5). These new z~5 
results build on earlier work by these authors (Iwata et 
al. 2003) and are the deepest and widest-area LF re- 
sults presently available at z~5, as they are drawn from 
a total of 1290 arcmin^ and reach MIqq+0.95, or 0.7 
mag deeper than the z^5 LFs of Ouchi et al. (2004). 
I. Iwata and his collaborators have kindly provided us 
with their new, deep estimates of z~5 incompleteness- 
corrected LBG number densities which they have com- 
puted using the Veff approach similar to the one we 
used in the KDF (see Iwata et al. 2003 and I. Iwata 
in preparation for details). We then used these num- 
ber densities to calculate the Schechter function parame- 
ters and their uncertainties following the same procedure 
as we used for our lower-z KDF samples (see KDF II). 
We find Mi*7oo = -20.80^°;^°, log0* = -3.10;°]^, and 
a ~ — I.OIIqIq at z^5 from the Iwata et al. data while 
the corresponding z~5 UV luminosity densities are given 
in Tables Hand 13 

Several interesting trends are apparent in Fig.O First, 
it is clear that the luminosity density at z~4 is substan- 
tially lower than that at ^^--^3. This is a direct conse- 
quence of the shallower faint-end of the luminosity func- 
tion that we find at 2;~4. This drop in luminosity den- 
sity at z^A is in contrast to the z^4 result of Steidel et 
al. (1999; small star symbols), but this is because these 
authors assumed a steeper, non-evolving a — —1.6 faint- 
end slope at this redshift. Our KDF results, instead, 
show an evolutionary trend whereby the luminosity den- 
sity increases with time from z^4 to z~3. As we show 
in more detail in § 13.31 this trend is due to the evolu- 
tion of the population of suh-L}^^Q galaxies. It is not 
clear whether this evolutionary trend continues to Z'^2, 



Thompson 

or whether there is a plateau in the luminosity density at 
0~3-l. Our data favor a plateau, but a peak at z^2 can- 
not be fully ruled out due to p otential systematic prob- 
lems at that redshift (see 12. 21 and Sawicki & Thompson 
2006 for details). 

It also appears plausible that the evolution we see at 
2^4— >3 in the KDF data may have started at earlier 
times: if we consider the Iwata et al. (2006) z~5 result 
then there appears to be a monotonic increase in the 
luminosity density produced by galaxies brighter than 
O.li^gQ that extends from very early times, at least 
z~5, until at least z~3. 

The GALEX points illustrate the now well-known drop 
in luminosity density from z~l to z~0. Together with 
our KDF results at z^A, 3, and 2, and the Iwata et al. re- 
sults at Z'^^^ the following trend is clear: the luminosity 
density, at least in galaxies brighter than O.IL^^q, in- 
creases steadily from the earliest redshifts until at least 
z~3; it then culminates somewhere between ^^--^3 and 
z~l as either a peak or a broad plateau; and it then drops 
precipitously from 2;~1 to 2^0. This picture is consis- 
tent with that first presented by Sawicki et al. (1997), 
who found a broad peak in the luminosity density at 
z^2.5 and an onset of cosmic star formation at z>3.5'^. 

3.2. The total luminosity density 

Our analysis thus far (summarized in Fig. |2Jl, has 
concerned itself with luminosity density due to galax- 
ies brighter than O.lL^^g. only. This is an approach 
common in the literature, but one that does not present 
a complete census of the cosmic luminosity density. We 
thus now turn to the question of how our results change if 
we incorporate the contribution of extremely faint galax- 
ies, L<L*j^gQ — galaxies that are too faint to be seen in 
any of the surveys to date at high redshift, including 
ours. 

The total luminosity density is obtained by integrating 
the luminosity function over all luminosities from to 
oo (Eq. This calculation involves an extrapolation 
to magnitudes fainter than our limiting magnitudes, but 
the depth of our KDF data allows this to be a robust 
extrapolation: (1) because of their large area and depth 
the KDFs have constrained the LFs' faint-end slopes suf- 
ficiently well, and, (2) the depth of our data is such that 
the extrapolation involves only ~25% of the total light 
(Fig.^l. Similar remarks apply to the surveys of Iwata 
et al. (2006) at z--5 and GALEX at z<l. Surveys with 
smaller areas and/or shallower limiting magnitudes must 
make much bigger and more uncertain corrections, but 
this is not the case for the KDF. 

The total luminosity densities are shown in Fig.^jwith 
open symbols and are compared with L>0.1L'^gQ re- 
sults of § 13.11 that are represented by filled points. The 
difference between the L>0.1L2bg ^^"^ total results 
is very small at high redshift, z>2. This is so because 
the faint-end slopes measured by us in the KDF (and by 
Iwata et al. 2006 at z~5) are relatively shallow, resulting 
in relatively little additional light in galaxies fainter than 
O.lL^gg. However, although the effect is small here, the 
importance of reliably measuring the faint end slope of 

^ We note that Madau et al. 1996 presented only lower limits on 
the high-redshift luminosity densities and that these values were 
corrected for incompleteness only later in Madau et al. 1998. 
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the LF cannot be understated: had the faint-end slopes 
turned out to be significantly steeper, much stronger dif- 
ferences would have ensued. As it is, most of the light at 
high redshift resides in galaxies with L=0.1-lL*]^gQ (as 
we show in S underscoring the importance of accu- 
rately measuring the LF's faint end slope at luminosi- 
ties fainter than i^BG- lower redshifts, z<l, there 
is a fairly significant difference between the total and 
the L>0.1L*j^gQ values. However, this is not because the 
low-z faint-end slope is very steep (it isn't), but primarily 
because at z<l L* is a magnitude or more fainter than 
our fixed high-z L} p,^. This difference ensures that the 
O.lL^gQ limit of S 13.11 does not encompass most of the 
low-z light. 

At z<l the total UV luminosity becomes dominated 
by very faint galaxies, galaxies with luminosities below 
O.IL^BG (^^^ Fig.E)). It is curious to note that the total 
luminosity density at z~0.5 is similar to that at 
which is not the case when the comparison is restricted 
to galaxies brighter than OAL^gQ only. Nevertheless, 
while the decline from z~l to z~0 is less precipitous 
when one considers the total luminosity density rather 
than just that due to galaxies brighter than O.IL^bG' 
the broad qualitative picture is the same in both cases: 
a gradual build-up in the luminosity density from very 
high redshift which is followed by a peak or a plateau 
between and z^^l which in turn is followed by a 
decline to z^O. 

This broad qualitative picture is also only mildly af- 
fected when one corrects for the effects of interstellar 
dust. Here, at each redshift we have applied a range 
of plausible dust corrections and the resulting dust- 
corrected star formation rate densities are shown as a 
gray band in Fig.^ We corrected the z>2 values by fac- 
tors of 5-15 as suggested by various studies of extinction 
in Lyman Break Galaxies (e.g., Sawicki et al. 1998; Pa- 
povich et al. 2001; Shapley et al. 2001, 2005; Vijh, Witt, 
& Gordon 2003; Reddy & Steidel 2004; Iwata et al. 2005) 
and the z<2 GALEX values by 2-5 (Schiminovich et al. 
2005; see also, e.g., Tresse et al. 2002; Perez-Gonzalez et 
al. 2003). After applying dust corrections we can now 
properly regard Fig. ^ as showing us the cosmic star for- 
mation rate history of UV-bright galaxies. 

Extremely dust-obscured galaxies such as high-z 
SCUBA sources (e.g.. Small et al. 1997; Barger et al. 
1998; Eales et al. 1999; Chapman et al. 2005) are not 
included in the dust-corrected star formation rate densi- 
ties shown in Fig. 01 Most such galaxies have UV lumi- 
nosities that heavily underrepresent their star formation 
rates and thus such galaxies are not properly included if 
our census is to be regarded as a measurement of star for- 
mation activity in the Universe. The number density of 
such heavily obscured sources at low redshift is too low to 
make a large contribution to the star formation rate den- 
sity. However, at high redshift, z~2, their numbers are 
sufficiently large that their contribution is significant and 
including them in our census would raise the star forma- 
tion rate density shown in Fig.^ The increase would be 
highest at z^2 where the redshift distribution of SCUBA 
sources appears to peak (Chapman et al. 2005). The cor- 
rection to the star formation rate density at that epoch 
may be up to a factor of two, i.e., 0.3 in dex (Chapman 
et al. 2005; Reddy et al. 2005). 

Taken together, however, the incorporation of dust cor- 



rections — be they only dust corrections to the UV- 
selected galaxies or also the inclusion of the contribution 
of heavily-obscured sources — does not qualitatively af- 
fect the picture shown in Figs.|21and0J Dust or no dust, 
the star formation rate density of the universe appar- 
ently began rising at some high redshift z>4, continued 
to rise until at least z~3, climaxed in a peak or a plateau 
between z~3 and z~l, and then dropped to z~0. 

While the decline from z~l to the present epoch has 
been well established for a long time (e.g., Lilly et al. 
1996; Shiminovich et al. 2005), the behavior of the cos- 
mic star formation rate denstity at z>3 was less clear: 
the Steidel et al. (1999) picture of a plateau extending 
from z~3 to z^A and possibly beyond (e.g., Giavalisco et 
al. 2004b; Ouchi et al. 2004) was for a long time the gold 
standard in the field. However, the Steidel et al. (1999) 
no-evolution result was based on the measurement of a 
non-evolving LF's bright end and the assumption of a 
non-evolving faint-end. Our KDF results clearly show 
that, in contrast to these assumptions, the luminosity 
function evolves and this evolution is reflected in the lu- 
minosity density which evolves with time at high redshift. 

3.3. Dependence on galaxy luminosity 

In the previous Section we have shown that the UV lu- 
minosity density exhibits a rise from early epochs to z~3, 
a peak or plateau between Z'--^3 and z~l, and then a de- 
cline to Z'-^^O. We now ask the question: which galaxies 
are responsible for this picture? To address this question 
we split our analysis into three luminosity intervals, com- 
puting the luminosity density in galaxies brighter than 

-^LBG' galaxies between L*j^bg ^'^'^ ^-^^lbg^ ^^'^ those 
fainter than O.IL^^q. 

The luminosity density due to just those objects with 
luminosities in the range L^ax > L > L„iin is given by 

PhiLmax > L> L,mn) = / L4>{L)dL (5) 

= </.*L*[7(a + 2,L„„„/L*) 
-"f{a + 2,L„iax/L*)\. 

We make the boundaries of our luminosity intervals at 
^LBG ^^'^ ^-^^LBG- l^ote that to compute the last 
(faintest) of those three intervals we must extrapolate 
the faint-end slope of the LF at all redshifts, but the 
data we consider here (the KDF, GALEX, and the z~5 
Subaru data of Iwata et al. 2006) all go sufficiently deep 
to measure the faint-end slopes well, thus making these 
extrapolations acceptable. 

In Fig. |S1 we plot the redshift evolution of the 
luminosity-dependent luminosity densities (colored sym- 
bols and lines) compared with the total luminosity den- 
sity from Fig. ^ in 13.21 fopen circles and black line). At 
high redshift, z>2, the luminosity density is dominated 
by galaxies in the narrow luminosity range {0.1~l)L^gQ: 
galaxies that are brighter than L*j^gQ or those that are 
fainter than O.IXJ^^q contribute relatively little light to 
the total at these redshifts! In contrast there is little evo- 
lution in the bright end of the galaxy population between 
z~2, 3, and 4 (and 5) and this results in a nearly con- 
stant luminosity density of bright, super-L^^^. galaxies 
over z~2-5. To date, most of the focus at high redshift 
has been on galaxies around L*j^gQ- But it is the sub- 
L* galaxies that, as we have illustrated here, dominate 
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the total luminosity density at z>Jl. This dominant sub- 
population is only now for the first time starting to be 
robustly probed by our KDF survey and by follow-up 
studies that are now either underway or being planned. 

We conclude this section by summarizing our main 
points. At high redshift, z>2, there is a rise in the to- 
tal luminosity density with time. This rise is not due to 
changes in the population of bright galaxies, L>L}^^q, 
but is instead directly related to the rise in the number 
density of faint, sub-L^^g, galaxies and the steepening of 
the faint end slope of the luminosity function that reflects 
it. 

4. DISCUSSION 

4.1. Comparison with other work 

In this section we compare our results with those of 
other surveys. Figure |3I shows results from several other 
programs that recently studied the cosmic UV luminos- 
ity density at high redshift. We consider here a number 
of recent surveys, but — with the exception of the work 
of Steidel et al. (1999) — we avoid results that are based 
on the very small area HDFs (e.g., Madau et al. 1996; 
1998; Sawicki et al. 1997; Connolly et al. 1997; Arnouts 
et al. 2005) because such results have to be held suspect 
given the small field of view of the HDFs and the result- 
ing susceptibility to cosmic variance. We have, however, 
included some z~6 results based on the small Hubble Ul- 
tra Deep Field since there are at present no other data 
that probe significantly below L*^qq at that redshift. In 
all the surveys that we do consider we have converted 
the published results to a consistent limiting depth of 
L=0.\L*^gQ and have adjusted them to rest- frame 1700A 
values as necessary. The superiority of the KDF results 
vis-a-vis these other surveys lies in the combination of 
the KDF's superb depth, large field of view, multiple 
sightlines, well-understood sample selection, and an ex- 
tensive range of robustness-gauging tests that we carried 
out as part of our LF measurement (Sawicki & Thomp- 
son 2006). 

z^2 and z^3 

At z~2 and z~3 our KDF results are in very good 
agreement with other available measurements, namely 
the photometric redshift results of Gabasch et al. (2004b) 
and the LEG work of Steidel et al. (1999). 

The agreement with Steidel et al. (1999) z~3 results 
is — at one level — not surprising because the KDF 
uses the very same filter set and sample selection as the 
ground-based Steidel et al. surveys. Indeed — our z^i 
and z^A luminosity functions incorporate Steidel et al. 
(1999) data at the bright end (see Sawicki & Thompson 
2006). On the other hand, however, Steidel et al. (1999) 
used the Hubble Deep Field to constrain the shape of 
their luminosity function below M£^q-|-1, thereby incur- 
ring the uncertainty of differential selection between their 
bright and faint galaxies, in addition to potential cosmic 
variance issues due to the HDF's small area. Although 
there is a small difference between our and their faint-end 
LF slopes at z~3 (our a=— 1.43 vs. their a=— 1.6), there 
is quite good agreement between the luminosity densities 

at i<0.1L2sG- 

Our KDF results are also in excellent agreement with 
the z-2 and z-3 FORS Deep Field (FDF) and Great 



Observatories Origins Deep Survey South (GOODS-S) 
resuhs of Gabasch et al. (2004b). The FDF is a sin- 
gle relatively small (~ 40 arcmin ) deep field with sam- 
ple selection (full-blown photometric redshifts) that is 
rather different than ours. Their GOODS-S data uses 
similar sample selection and area (~ 50 arcmin^) as the 
FDF but is a magnitude shallower than their FDF data, 
reaching only ^M^^q+I at z^'i — depths comparable 
to those reached from the ground by Steidel et al. (1999). 
Despite the differences in sample selection and the rel- 
atively small area of their significantly deep data (the 
FDF), there is very good agreement between the lumi- 
nosity densities they find and those that we find in the 
KDF. 

OveraU, aU three results (KDF, Gabasch et al. 2004b, 
and Steidel et al. 1999) are highly consistent with each 
other, albeit ours is the most robust of the three given 
the KDF's large area and multiple sightlines (important 
for mitigating the effects of cosmic variance) combined 
with faint limiting magnitude (vital for constraining the 
faint end of the galaxy population), and extensive LF 
robustness tests. 

At Z'^A, we compare our results with Subaru Deep Sur- 
vey (SDS) work of Ouchi et al. (2004), Giavalisco et al. 
(2004b) results based on the Great Observatories Ori- 
gins Deep Survey (GOODS; Giavahsco et al. 2004a) HST 
data, as well as the aforementioned work of Gabasch et 
al. 2004b and that of Steidel et al. (1999). At z-4, the 
GOODS-S component of the Gabasch et al. (2004b) anal- 
ysis reaches only ~Af£^g,-|-0.5, so here we consider their 
work to be based primarily on the single-pointing FORS 
Deep Field. All the surveys considered use variants of 
the Lyman Break color-color selection technique, with 
the exception of the FDF, which uses full photometric 
redshifts. 

These four surveys (Ouchi et al. 2004; Gabasch et al. 
2004b; Giavahsco et al. 2004; and Steidel et al. 1999) 
do not present a consistent picture, giving a range of 
luminosity densities that span a factor of ^--^2.5. Let us 
examine the differences in more detail. 

Gabasch et al. (2004b) do not give FDF luminosity 
density measurements at z~4 but present values at z=3.5 
and 4.5 (see Fig. OJ. However, these values bracket our 
z^A KDF result and a straightforward interpolation be- 
tween their FDF z=3.5 and z=A.5 points is in very good 
agreement with our data. This agreement is to be ex- 
pected given the good agreement between our LFs and 
theirs (see Sawicki & Thompson 2006), but note that our 
KDF result is much more robust for reasons outlined in 

§Eini 

The results of Steidel et al. (1999) at first appear to 
be somewhat higher than ours, but this is because of 
their assumption that the faint end of the LF does not 
evolve from z~3. With a shallower faint-end a at z~4, 
as found in the KDF, their results would be in agreement 
with ours. 

Giavalisco et al. (2004b) give two measurements from 
the GOODS data, obtained using two different tech- 
niques. One measurement (from their technique) 
paints a picture of a luminosity density that docs not 
differ from that at z~3, while the other (from their Veff 
approach) yields an evolving luminosity density, in agree- 
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ment with our results. Giavalisco et al. prefer their 
method and results, but our KDF data are in much bet- 
ter agreement with their lower, 14// point. 

The SDS work of Ouchi et al. (2004) gives the highest 
luminosity density value but this value may well be prej- 
udiced by potential problems with their measurement of 
the LP's faint end: at their data probes one mag- 
nitude fainter than M^^q and they find an extremely 
steep faint-end slope of a=— 2.2 that gives an infinite to- 
tal luminosity density (Eq. A strong limitation of 
their work is that their color selection is not well cali- 
brated spectroscopically except for a handful of redshifts 
and must instead rely heavily on models of galaxy col- 
ors. It is not clear how strong are the systematic effects 
inherent in this modeling. 

Overall, our results are in very good agreement with 
those from the PDF and with the GOODS K// ap- 
proach, while we find substantially less luminosity than 
does the SDS and the GOODS approach. We feel 
that, given the robustness of our sample selection and 
LP measurement, our results are the most reliable of all 
those currently available at z^A. 

4.2. Earlier epochs 

The literature now contains several measurements of 
the luminosity function and cosmic luminosity density 
at z~5. Here we discuss those by Iwata et al (2006) as 
well as those based on GOODS (Giavalisco et al. 2004b), 
SDS (Ouchi et al. 2004), and PDF (Gabasch et al. 2004a, 
2004b; note that the PDF results are at Z--4.5). With 
the exception of the Giavalisco et al. GOODS results, 
all these surveys present both luminosity densities and 
luminosity functions; of these, the Iwata et al. (2006) 
Subaru work has the best combination of depth (0.7 mag 
deeper than the SDS) and area (an order of magnitude 
larger than the PDF). 

The Iwata et al. (2006), PDF, and the GOODS K// 
results all indicate that the luminosity density at z~5 is 
lower than that at z~4 in the KDF. In both the Subaru 
work of Iwata et al. (2006) and in the PDF (Gabasch 
et al. 2004a) the faint-end slope of the LP is quite shal- 
low at these redshifts, in apparent continuation of the 
LP evolutionary trend that we see in the KDF. It is this 
shallow faint-end slope that results in the low luminosity 
density that these surveys measure at z^5. On the other 
hand, the GOODS results as well as those from the 
SDS would suggest a higher luminosity density — one 
that is almost unchanged from z~3. In the SDS work 
this is caused by a very steep faint-end slope that they 
assume from their lower-z results (their observations at 
z~5 reach only to LI^q+0.25). The z~3^4 trend seen 
in the KDF tends to favor a continuing drop in the lumi- 
nosity density to higher redshifts, in line with the Iwata 
et al. (2006), PDF, and GOODS-K// values. 

zr^6 

The situation at z~6 is tenuous despite much recent 
effort. Some recent results (Bouwens et al. 2006; the 
GOODS results of Giavalisco et al. 2004b) suggest 
a z^6 luminosity density nearly as high as that at z~3. 
Others, (Bunker et al. 2004; Ve// results of Giavalisco 
et al. 2004b) give values that are significantly lower (see 



Fig.EJ. 

The situation is tenuous for several reasons. First, 
at z~6 sample selection is generally done using just 
one color, i—z, which makes it difficult to estimate red- 
shift distributions, and hence survey volumes, and hence 
source number densities. The absence of a third, redder 
band also makes source luminosities very uncertain be- 
cause at z~6 the z-band fluxes that are invariably the 
starting point for computing luminosities are attenuated 
by intergalactic hydrogen line blanketing. While correc- 
tions for this attenuation are usually applied, such cor- 
rections are necessarily statistical because they rely on 
an average rather than exact distribution of intervening 
gas, and so are subject to systematic biases that are dif- 
ficult to quantify. Moreover, these corrections are typi- 
cally based on the models of Madau (1995); these models 
are based on lower-redshift data, and it is unclear how 
well they work at z~6. It is thus not clear how well 
the fluxes and hence luminosities of z^6 candidates are 
being estimated. 

Another limitation is that the all-important shape of 
the faint end of the z^6 luminosity function is highly 
uncertain. Assuming the Steidel et al. (1999) z~3 faint- 
end slope of a=— 1.6 is not a good approach in hght of 
the fact that the true slope is shallower even at z~3 and 
that it becomes progressively shallower with increasing 
redshift (Sawicki & Thompson 2006; see also Gabasch et 
al. 2004a). On the other hand, empirically constraining 
the number density of sub-L^^g. galaxies from the data 
is highly problematic since the only field sufficiently deep 
for this task is the Hubble Ultra Deep Field (UDF). The 
UDP's small area (12 arcmin^) makes it subject to poten- 
tial cosmic variance effects only slightly less severe than 
those affecting the HDPs. Bouwens et al. (2006) argue 
that cosmic variance effects on UDF scales are not large 
at Z'--^6 (they suggest a 35% RMS effect), but this seems 
questionable: e.g., the difference in comoving luminos- 
ity density at redshifts Z'^S and Z'^A between the two 
~5 arcmin^ Hubble Deep Fields is a factor of two (e.g., 
Ferguson, Dickinson, & Williams 2000) and it is reason- 
able to suppose that similarly large fluctuations can be 
present on the just slightly larger scale of the UDF. 

In summary, present measurements of the z~6 UV lu- 
minosity density are highly uncertain. Given that our 
KDF data show that luminosity density is dropping al- 
ready from z~3 to z~4, it seems unlikely that it would 
be as high at z~6 as it is at z~3. A simple extrapolation 
of the trend we see between z~3 and z^-^A in the KDF 
— a trend that seems to continue to z~5 in the work of 
Iwata et al. 2006 — suggests the cosmic UV luminosity 
density at z~6 is similar to that reported by Bunker et 
al. (2004). To make further progress on this issue will re- 
quire large, multi-field z~6 galaxy samples with J-band 
photometry that reaches significantly below L^bq- 

4.3. Keeping the Universe ionized 

Recently there has been much interest in the ques- 
tion of the objects that are responsible for reionizing 
the Universe at high redshift (e.g., Lehnert & Bremer 
2003; Bunker et al. 2004; Martin & Sawicki 2004; Sti- 
avelfi. Pall, & Panagia, 2004; Yan & Windhorst 2004). 
Several of these authors have concluded that the num- 
ber of luminous galaxies at z^6 is insufficient to provide 
enough photons to keep the Universe ionized and have 
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postulated that sub-L* galaxies that are largely too faint 
to be observed directly at these redshifts may make up 
the deficit. The KDF's highest redshift bin is at z~4, 
only 0.7 Gyr after the putative redshift of reionization 
at 2=6.5, and since our survey probes significantly be- 
low L*j^BGy ^ good position to comment on the 
issue. We will study here whether faint galaxies produce 
enough ionizing photons to maintain ionization at Z'^A 
and later, and ask what that tells us about keeping the 
Universe ionized at earlier epochs. 

Following much of the z^6 work (e.g., Ferguson, Dick- 
inson, & Papovich 2002; Bunker et al. 2004; Yan & Wind- 
horst 2004; Bouwens et al. 2006), our estimate will be 
based on the Madau, Haardt, & Rees (1999) radiative 
transfer calculation. As we mention later, this is not 
necessarily the best approach to adopt, but it is the ap- 
proach taken by most observers so for consistency we, 
too, will use it here. 

According to Madau, Haardt, & Rees (1999), the rate 
of production of ionizing photons needed to maintain ion- 
ization at a given redshift can be expressed as 



'1.58 • 10^1 



(6) 



/1.58-lO^S^ /l + zx3/17fc/i^„y 

I s-Mpc3 J^3°l~^J lo:o4rJ 



where ni/(0) is the average density of hydrogen, trec{z) 
is the average recombination time, C30 is the clumping 
factor of ionized hydrogen, f2f, is the baryon density, and 
/170 is Hubble's constant in units of 70 km s^^Mpc^^. It 
should be noted that Madau et al.'s 6*30=1, motivated 
by simulations by Gnedin & Ostriker (1997), is likely 
too high and more recent studies suggest reionization 
is somewhat easier to achieve (see, e.g., Miralda-Escude 
2003; Furlanetto & Oh 2005; Meiksin 2005). At rest- 
frame I700A, we are measuring photon energies that are 
too low to ionize hydrogen and so we need to relate our 
I700A luminosity densities to number densities of A=0- 
912A photons. Using the Starburst99 evolutionary spec- 
tral synthesis models (Leitherer et al. 1999) we find that 
the rate of production of ionizing photons. A/ion, is re- 
lated to the I7OOA luminosity ^1700 by 
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Equation [7| is valid for starbursts with age >10^yr, 
Salpeter (1955) stellar initial mass function (IMF) with 
1OO>M/M0>O.1, and is insensitive to metallicity. Then, 
combining EquationsEland[3we find that if the Universe 
is to be kept ionized by ongoing star formation, we need 
to have a 1700A luminosity density of no less than 



-^1700 (^) 



erg • s ^Hz ^Mpc 



9.63-10 
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I 
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(8) 

We have included in Eq.|Sla relative escape fraction term 
/Js^, which accounts for the fact that in practice, due 
to various radiative transfer effects, fewer ionizing pho- 
tons escape a galaxy than would be straightforwardly 
predicted by the number of escaping non-ionizing pho- 
tons. There is currently much uncertainty about the es- 
cape fraction of ionizing radiation, but it seems reason- 
able that /e7i<0.5 (Heckman et al. 2001; Steidel, Pettini, 



& Adelberger 2001; Fernandez-Soto, Lanzetta, & Chen 
2003; Inoue et al. 2005), and certainly ^7^<1 since the 
relative escape fraction cannot exceed unity. 

The amount of UV starlight needed to keep the Uni- 
verse ionized (as predicted by Equation IHJ is shown as a 
solid gray curve in Figs. |21 and 0] for hydrogen clumping 
factor value of C3o=l and a conservative escape fraction 
/esc=0-5- Lower escape fractions would push the curve 
upwards while lower clumping factors, C3o<l, would 
lower it. As has been pointed out by a number of au- 
thors, z^6 galaxies brighter than O.lLJ^^g do not pro- 
duce enough UV light to ionize the Universe (Fig. OJ. 
However, it would seem that there isn't enough starlight 
to keep the Universe ionized even at z~4, and this is so 
not only if we restrict ourselves to galaxies brighter than 
0.lL*j^QQ, but even if we integrate down the luminosity 
function all the way to L—Ol 

In the calculation presented here, an unevolving z~4 
galaxy population would clearly not be sufficient to main- 
tain ionization at Z'^Q, and yet our observations indicate 
that the galaxy population is evolving, with the UV lumi- 
nosity density dropping with increasing redshift. While 
at z~4 additional radiation from active galactic nuclei 
(AGN) might be enough to make up the deficit (but 
see Steidel et al. 2001), this is unlikely to be the case 
at higher redshifts as the number density of even low- 
luminosity AGN is decreasing with increasing lookback 
time (e.g., Willott et al. 2005). It would thus seem that 
the Universe cannot produce enough ionizing radiation 
either from stars or black holes to reionize its gaseous 
content at z~6. 

There are several ways out of this apparent paradox 
but ionizing radiation from very faint but otherwise "nor- 
mal" star-forming galaxies does not appear to be one 
of them. Some possibilities include a radically non- 
standard IMF or stars with very low, nearly primordial 
metallicities (e.g.. Bunker et al. 2004; Stiavelli, Fall, & 
Panagia 2004). A less radical alternative, however, is 
that the Madau et al. (1999) requirement for reioniza- 
tion that is assumed by most observers is overly strict. 
There is an ongoing debate about the number of pho- 
tons needed for reionization, a debate that is related 
to the uncertainty in the dumpiness of the intergalactic 
medium that controls the importance of recombinations 
(e.g., Madau et al. 1999; Miralda-Escude 2003; Furlan- 
etto & Oh, 2005; Meiksin 2005). It now seems highly 
plausible that the C3o=l assumed by Madau et al. (1999) 
and subsequently by many observers is significantly too 
high. If, for example, we follow the results of Meiksin 
(2005) instead of those of Madau et al. (1999) then the 
required L^^qq is a factor of ~5 lower than that shown 
in Figs, inland^ Under these less restrictive clumping 
assumptions massive stars in faint galaxies that follow 
the even relatively shallow faint-end LF slope we find in 
the KDF can produce more than enugh ionizing photons 
at z~4 and plausibly also at higher redshifts. In this 
case, there is no need to invoke a very numerous faint 
population and steep LF faint-end slope at 2:^6. 

4.4. The galaxies that dominate the UV luminosity 
density at high redshift 

The steep increase in the total luminosity density with 
time at high redshift is caused by the increase in the 
number density of sub-L^^^ galaxies and the associ- 
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ated steepening of the LF's faint-end slope. However, 
while the faint end of the luminosity function evolves, 
the bright end stays nearly constant — a fact that is 
reflected in the relative flatness of the luminosity den- 
sity due to L>L*]^gQ galaxies. This differential evolution 
suggests that individual galaxies themselves are evolving 
and that they are doing so differentially with luminosity. 

The evolutionary behavior of the differential luminos- 
ity density that we see over z~5— s-O (Fig. ISJ is reminis- 
cent of the galaxy downsizing picture (e.g., Cowie et al. 
1996, Iwata et al. 2006a). Here, the epoch of vigorous 
star formation in the most luminous galaxies — those 
with L>L2bg — ends first with a rapid decline in their 
activity starting at z~2. This is followed at later epochs 
by a drop in the activity due to intermediate-luminosity 
galaxies (L=0.1-li2gQ). This in turn is followed by the 
dominance of those galaxies in our faintest luminosity 
bin — LjO.lL^sG- However, unlike in the lower-redshift 
work of Cowie et al. (1996), here, at higher redshifts, the 
downsizing picture is so far only seen in luminosity and 
it remains to be seen if it is also tied to galaxy masses. 

The differential evolution we observe in the luminos- 
ity density and the luminosity function must reflect evo- 
lution of individual galaxies. However, at present the 
mechanism for this evolution remains unknown. In Saw- 
icki & Thompson (2006) we have suggested several possi- 
ble mechanisms related to changes in timescales of star- 
forming episodes, or evolution in the amount or distri- 
bution of their interestellar dust. It is also possible that 
we are witnessing two different evolutionary mechanisms 
that dominate the different ends of the luminosity func- 
tion and are competing around the characteristic LF lu- 
minosity L* . It is possible, for example, that on the 
two different sides of L* we are seeing galaxies that are 
governed by two different regimes of feedback — e.g., 
AGN feedback in more massive (and plausibly more lu- 
minous) galaxies, and star- formation-driven feedback in 
lower-mass (and lower luminosity) systems. 

We are now working to understand the nature of the 
evoltion we see by studying differences in the properties 
of individual galaxies as a function of epoch and luminos- 
ity. Clustering and spectral energy distribution studies 
will let us constrain the dark halo and stellar masses of 
galaxies as a function of luminosity, and thus will tell us 
whether the downsizing we see is related to galaxy masses 
as it seems to be at lower redshifts (Cowie et al. 1996). 
SED studies will also let us test whether other galaxy 
properties, such as their dust properties or timescales 
of their starbursting episodes vary with luminosity and 
epoch and thus whether they are responsible for the evo- 
lution we see. While observationally expensive, such dif- 
ferential comparisons of faint members of a population 
with their brighter brethren can lead to important in- 
sights (see, e.g., Sawicki et al. 2005) and will be key to 
understanding what drives the evolution of galaxies at 
high rcdshift. 

5. SUMMARY AND CONCLUSIONS 

Our KDF survey has allowed us to make the best mea- 
surement thus far of the UV luminosity density at z~4, 
3, and 2. This is the case because our work builds di- 
rectly on the well-tested photometric selection techniques 
of Steidel et al. (1999, 2003, 2004), because it probes very 
deep into the galaxy luminosity function, and because it 



does so using several spatially independent fields to con- 
trol for the effects of cosmic variance. The main results 
of our analysis of the UV luminosity density are: 

1. The picture of total UV luminosity density, i.e., lu- 
minosity density due to galaxies of all luminosities, 
is that of a rise from early epochs, z>4, followed 
by either a plateau or a broad peak at 2;~3-l, and 
then followed by a decline to z=0. This picture 
is consistent with that presented by Sawicki et al. 
(1997), although is at odds with some recent work 
that posits a constant luminosity density over Z'^G— 
1. 

2. The luminosity density at high redshift, ^^--^2-4, 
is dominated not by luminous galaxies, but by 
snb-L*j^gQ galaxies, where L*^^q is the characteris- 
tic 1700A luminosity that corresponds to M^^q— 
—21.0. Indeed, it is galaxies within a rather narrow 
luminosity range L=Q.\-\L*^gQ that dominate at 
these redshifts and it is the rapid evolution of these 
galaxies that drives the evolution in the total UV 
luminosity density. 

3. In contrast to 2:>2, at lower redshifts, z<l the to- 
tal UV luminosity density, as shown most recently 
by GALEX, is dominated by very faint galaxies 
with L<Q.lL*j^gQ. The L=0.1-1L2bg galaxies 
that dominated at z>2 are still important but no 
longer dominant. This reversal in the luminosity 
of the galaxies that dominate the total luminos- 
ity density is reminiscent of galaxy downsizing, al- 
though the effect we see so far is a downsizing in 
UV luminosity rather than galaxy mass. 

4. Several recent studies have suggested that at z^6 
there are not enough luminous galaxies to maintain 
ionization of intergalactic gas and proposed that 
a large population of faint galaxies is needed for 
that task. Within the context of the most widely- 
used models, we find that there seemingly aren't 
enough bright and faint star-forming galaxies even 
as late as 2^4 to maintain ionization, questioning 
the assumptions of these most widely-used models. 
Apparently the Universe must be easier to reionize 
than some recent studies have assumed. 

In this paper we have shown that the total UV lumi- 
nosity density of the Universe at high redshift is dom- 
inated by galaxies fainter than the L*^gQ galaxies that 
have thus far drawn most of the community's attention. 
These sub-L^^g, galaxies are, as we have also found in 
Sawicki & Thompson (2006), evolving rapidly as a pop- 
ulation, while the population of their far better stud- 
ied, more luminous cousins appears to remain unchanged 
with time. The question then becomes, what makes the 
high-z galaxies brighter than L'^gQ and those fainter 
than L*j^gQ so different. We will address this issue in the 
future by studying the clustering properties and stellar 
populations of these important systems. 

We thank the time allocation committee for a gener- 
ous time allocation that made this project possible and 
the staff of the W.M. Keck Observatory for their help 
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Fig. 1. — The amount of total light above a limiting magnitude Mn^in a Shechter function with faint end slope a. The colored lines 
represent our KDF results at z~3 (green) and (orange) and the ground-based z-^S and z~4 surveys of Steidel et al. (1999; blue). The 
solid parts of these lines show the directly observed galaxies (down to the respective limiting depths of these surveys) and the dashed lines 
illustrate extrapolations below the completeness limits. The Steidel et al. (1999) ground-based work detects only ~l/3 of the UV light at 
z~3 (and even less at z~4), while the KDF captures ~75% of the UV light at these redshifts. 
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Fig. 2. — The luminosity functions of UV-selected star-forming galaxies at 2:~4, 3, and 2.2 derived by Sawicki & Thompson (2006). The 
z~4 and z~3 LFs are very robust against systematic effects; at z~2.2, the number density of galaxies is not less than that shown and if it 
is higher than shown, then it is so by no more than a factor of ~2. 
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The curve labelled i^TOO shows the amount of 1700A light expected — on the basis of the calculation in ? I4.3I — to be needed if starforming 
galaxies are to maintain ionization of the intergalactic gas. 
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Fig. 4. — The total luminosity density and dust corrections. Filled symbols show the luminosity density due to galaxies brighter than 
L'j^gQ, while open symbols show the total (all L) luminosity density. The gray band shows the allowable range after the light from the 
total (all L) UV-bright galaxy population has been corrected for dust obscuration as discussed in § 13.21 
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Table 1. Rest-frame 1700 A UV luminosity densities* 



z 




total 






L<0.1L*ac; 


survey 


0.05 


25, 


cc+0.08 
■'-"J—O 08 


18.56l|;« 


24.74 


+0.14 
-0.14 


25.49 tT& 


GALEX 


0.3 


25, 


04+0:17 


9n 7-3+4.90 


25.24 


+0.25 
-0.25 


25.72 


GALEX 


0.5 


25, 


■^"-0.40 


23.54l'^;g 


25.48 


+ 0,29 
-0.29 


25.78 


GALEX 


0.7 


26, 


14+0.56 


24.19l°jl 


25.69 


+0.32 
-0.32 


25 94 +0-66 
zo.at _o.66 


GALEX 


1.0 


26, 


19+0.73 


24.45t0;|° 


25.69 


+0.33 
-0.33 


25 91 +°-^'' 


GALEX 


1.7 


26, 


7c:+0.02 
■ '^^-0.02 


25.5910:11 


26.63 


+0.02 
-0.02 


96 03 +0 05 
_o.05 


KDF 


2.2 


26, 


40+0.03 
■^■^-0,03 


25 fi6+"-°^ 


26.21 


+0.02 
-0.02 


25 76 


KDF 


3.0 


26, 


49+0.03 
■ -^^-0.03 


25 64+"-"^ 


26.15 


+0,02 
-0.02 


25 91 +0-10 


KDF 


4.0 


26, 


Qg+0.07 
■"°-0.07 


25 47+0-05 


25.81 


+0.05 
-0.05 


25.34 t0;2| 


KDF 


4.8 


25, 


og+0.09 


9c; 0-3+0.08 


25.63 


+0.09 
-0.09 


24.93 


Iwata et al. (2006)'' 



*In logarithmic units of erg/s/Hz/Mpc"^ 

''Computed by us from data kindly provided by these authors in advance of publication 



Table 2. UV luminosity densities in galaxies brighter than O.IL* 



Survey 


2 


log 


;(pc/v"-) 
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0.05 


24. 


74+0,14 

'^-0,14 




0.3 
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■^^-0,25 




0.5 
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.0+0.29 

■'*°-0,29 




0.7 


25, 


7n+0,34 
'^-0,34 




1.0 
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yi +0,35 
'-'-0,35 


KDF 


2.2 
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-,9+0,02 
■^^-0,02 




3.0 


26, 


9R+0.01 
^"-0.01 




4.0 


25, 


Q7+0.03 

^' -0,03 


Iwata et al. (2006)'' 


4.8 


25, 


01 +0,07 
"-^-0,07 


Steidel et al. (1999) 


3 


26. 


-in+O-Ofi 
■■="^-0,06 




4 


26. 


99+0,10 
^^-0,10 


Giavalisco et al. (2004) method 


3.78 


26. 


-,q + 0.05 

•^"^-O.OS 




4.92 
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5.74 
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10+0. 24 
'°-0,19 


Giavalisco et al. (2004) 14// method 


3.78 


26. 


19+0.05 
^^-0,05 




4.92 


25. 


ot;+0.08 
'''^-COS 




5.74 


24. 


Q9+0.24 
^^-0.19 


Bunker et al. (2004) 


6 
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cc+o.oe 

"^■^-0,06 


Bouwens et al. (2006) 


6 


26, 


in+o,05 
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Ouchi et al. (2004) 
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-30+0,04 
■J"-0,04 




4.7 
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Gabasch et al. (2004a, 2004b) 
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3.5 


26, 
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4.5 
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''Rest-1700A luminosity density in units of erg/s/Hz/Mpc^ 

''Computed by us from data kindly provided by these authors in advance of publication 



